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Abstract 

Energetic  implosions,  using  two  or  three  load  wires 
to  create  a  focused  axial  stagnation  of  dense  wire  cores 
amidst  the  assembled  precursor  plasma,  are  examined 
with  respect  to  the  trade  between  the  implosion  mass 
lost  to  precursor  ablation  and  the  mass  or  kinetic  energy 
available  at  stagnation.  The  calculated  kinetic  energy  at 
stagnation  serves  as  measure  for  the  output  x-radiation. 


L  INTRODUCTION 

The  use  of  a  single  wire  load  on  large  machines  like 
"Z"  is  usually  avoided  due  an  excessive  initial  inductance 
and  subsequent  high  voltages  on  upstream  components. 
There  is  however  a  middle  ground  where  a  few  larger 
wires  can  start  well  off  axis  and  present  much  less  initial 
inductance.  Such  loads  will  provide  the  lower  inductance 
at  the  expense  of  more  precursor  plasma  involvement. 

In  contrast  to  closed  arrays  with  hundreds  of  fine 
wires,  this  relatively  unexplored  path  to  energetic  im¬ 
plosions  would  use  two  or  perhaps  three  load  wires  of 
appropriately  heavier  mass  and  aims  to  create  a  focused 
axial  collision  of  dense  wire  cores  amidst  the  assembled 
precursor  plasma.  The  precursor,  not  confined  by  many 
wires,  would  presumably  not  soften  this  collision. 

Hence  we  must  examine  the  trade  between  the  implo¬ 
sion  mass  lost  to  precursor  ablation  and  the  mass  avail¬ 
able  to  deliver  the  wire  core’s  kinetic  energy. 

A  clear  consequence  of  this  load  choice  is  a  quite  open 
geometry  for  which  the  transition  to  a  highly  conducting 
annular  MHD  plasma  is  neither  an  early  nor  a  necessarily 
dominant  feature  of  the  electrodynamics.  Generalizing 
slightly  beyond  MHD,  we  use  a  Lorentz  gauge  direct 
field  solver  to  treat  the  TEM  to  TM  mode  set  transition 
in  the  pinch  region.  The  problem  is  formulated  using 
a  scalar  potential  <f>  and  an  axial  vector  potential  Az 
as  unknowns,  with  the  radial  vector  potential  Ar  being 
determined  by  the  gauge  constraint. 


In  contrast  to  earlier  wire  dynamic  model  (WDM) 
formulations  with  inductance  matrix  elements  good  only 
in  the  thin  wire  limit  [1],  the  present  work  will  make  use 
of  a  new  analytic  Green’s  function  for  Az  that  accounts 
for  proximity  effects  among  the  wire  cores  and  the  return 
current  structure. 

The  Green’s  function  provides  a  first  estimate  for  the 
axial  vector  potential  solution  Az(x,  j)  at  the  pinch  mid¬ 
plane  and  the  net  load  inductance  arising  among  the 
heavy  wire  filaments.  That  solution  is  further  refined  by 
a  current  density  source  term  representing  the  extended 
coronal  plasma  conductivity  and  flow. 

The  use  of  empirical  mass  source  terms  [2],  for  estimat¬ 
ing  the  precursor  mass  ejection  rate  and  ablation  velocity, 
is  then  introduced  to  complete  the  picture. 


II.  MATHEMATICAL  FORMULATION 

There  are  two  required  innovations  to  get  the  results  we 
want.  First,  we  must  go  beyond  the  thin  wire  induc¬ 
tance  limit,  which  is  fine  for  wire  core  dynamics  and 
current  sharing,  but  starts  to  break  down  when  the  tenu¬ 
ous  corona  plasma  jackets  the  wires  and  entrains  current. 
Moreover,  in  the  larger  wires  used  for  the  loads  we  con¬ 
sider,  proximity  effects  become  noticeable  as  the  wires 
approach  one  another.  Hence  a  means  of  describing  the 
extended  conduction  medium  of  the  corona  while  main¬ 
taining  contact  with  the  lowest  order  inductive  picture 
enforced  by  the  wire  cores  is  key  to  a  versatile  picture. 

Second,  we  require  a  fluid  description  that  admits  a 
large  adaptability  in  scale  lengths  in  order  to  track  the 
history  of  coronal  plasma  elements  that  drift  radially  to 
form  the  axial  precursor.  A  smooth  generalization  of  the 
discrete  wire  filaments  into  fluid  particles  allows  a  large 
dynamic  range  in  mesh  size  for  the  field  solver.  The 
new  fluid  particle  methods  used  here  adapt  easily  to  such 
grids  because  the  projection  of  the  fluid  particle  variables 
like  density,  velocity,  and  energy  onto  the  Eulerian  grid 
is  virtually  exact. 
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A.  Energy  Transfer  in  a  Circuit  Model  B.  Field  Methods  for  Open  Geometries 


For  the  original  WDM  the  direct  interaction  among  an 
ensemble  of  current  elements  is  equivalent  to  a  potential 
energy  defined  by  the  array  inductance.  Within  the  en¬ 
semble  are  symmetry  groups  of  multiplicity  5,  e.g. 

/=CV',*=^’  (1) 

where  a-  is  any  5-wise  invariant  current  fraction  con¬ 
tained  in  any  particular  wire  path  represented  by  a  series 
resistance  and  inductance. 

By  construction  then  the  sum  of  all  a-  is  one.  Since 
the  wire  paths  are  in  parallel,  the  voltage  seen  by  each 
path  a-  is  equal  in  the  absence  of  (small)  wave  transit 
time  effects.  If  we  examine  the  voltage  Ve(t)  impressed 
at  the  entrance  or  feed  of  a  wire  array  cage,  then  from  the 
following  figure,  the  resolution  of  the  set  a-  and  the  total 
current  I  admitted  to  the  array  are  seen  to  be  two  distinct 
problems.  From  the  view  of  the  external  generator  Vg(t), 
the  array  is  a  two  terminal  device  and  the  energy  into  it 
must  be  the  same  for  all  sets  a-  that  show  the  same  net 


Figure  1.  WDM  circuit  elements  can  decompose  into  symmetry  groups, 
with  each  element  representing  a  resistance  and  inductance. 


Maintaining  equal  voltage  at  each  circuit  element,  the 
equations  fixing  the  current  fractions  a-  for  any  given 
Ve(t)=Vg(t)  —  ZgenI  and  I  can  be  written  as, 

Ve{t)  =  Ri(t)ail+{jiLirajI^  ,  (2) 

and  then  integrated  over  a  time  interval  <5  =  W7! T,  t0 
eliminate  the  time  derivative.  Denote  by  A¥e  the  change 
in  flux  over  this  interval,  the  matrix  relation 

~fui  SRjltj  CCj  •  /  •  otj  (3) 

emerges,  with  I ^  the  diagonal  identity  matrix,  u-  a  unit 
vector  over  the  local  group. 

In  the  limit  of  vanishing  inductance,  or  slow 
timescales  (5  Li/Ri ),  it  is  easy  to  see  that  the 

a-  =  Z, , //?■(?)  and  thus  they  clearly  add  up  to  unity. 


Loads  with  fewer  wires  using  larger  initial  wire  diam¬ 
eters  and  load  radii  will  present  a  lower  initial  inductance 
at  the  expense  of  more  precursor  plasma  involvement.  In 
contrast  to  closed  arrays  with  hundreds  of  fine  wires,  this 
very  open  field  geometry  is  positioned  on  the  edge  of 
validity  for  a  conventional  inductance  and  MHD  picture 
due  to  the  large  dynamic  range  in  scale  lengths  and  mag¬ 
netic  Reynolds  number.  The  limit  of  spatially  constant 
voltage  throughout  the  corona  is  probably  pretty  good 
but  by  no  means  certain. 


7.  Mode  Transitions  from  TEM  to  TM 


With  ZQ  —  \f  £0f  ll0  —  367.70,  an  inlet  TEM  boundary 
condition  on  the  scalar  potential  0(r,  z,  t) 


Vin=ZJc 


(4) 


determines  all  the  fields  near  the  inlet,  where  A  — »  0. 

Elsewhere,  the  time  integral  of  the  scalar  potential 
forms  a  useful  generalization  of  the  familiar  inductive 
“flux  function”,  viz.  let  xF(r,z,?)  =  and 


>P(r>) ->P(r<)  =  /"  A  r{zc,rvt)dry,  (5) 

Jr< 

'p(z>)-'p(z<)  =  [  y  Az{zvrc,t)dzl  ,  (6) 

Jz< 

on  radial  zc  and  axial  rc  conductor  boundaries  to  ensure 
that  the  tangential  component  of  electric  field  vanishes. 


2.  Utility  of  the  Lorentz  Gauge 

The  Lorentz  gauge  condition  choice 

r~ldr(rAr)  +  dzAz  +  c~2dt^>  =  0  ,  (7) 

allows  the  dynamics  to  concentrate  on  the  solution  of 
A Z(z>  rf  )  only.  One  may  in  fact  examine  only  two  cou¬ 
pled  wave  equations  for  the  needed  fields. 

V20-c“2^20  =  0,  (8) 

V2Az-c-2c>2  Az=-iUz,  (9) 

with  the  axial  current  density  near  the  wires  given  by 

Jz  =  a±[Ez  +  VrxBe(A)]  .  (10) 

For  the  inhomogeneous  wave  equation  involving 
A z(x,y),  a  Green’s  function  can  be  developed  for  the 
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problem  at  the  diode  mid-plane  which  captures  the 
proper  boundary  conditions  on  the  outer  return  current 
surface  and  on  the  finite  size  conducting  wire  cores. 
The  function  can  be  built  up  piecewise  from  image 
currents  generated  by  a  source  cylinder  (a),  any  number 
of  floating  cylinders  b,b’,. . .  and  the  return  cylinder  (c), 

&{x,y\s  =  a,b,b' . . .)  =  ga{x,y)  +  gb{x,y)  + 

g»{x,y)  +  gxb(x,y)  +  gxb,(x,y)  +  . . .  ,  (11) 

and  over  the  symmetry  groups  A z(x,y)  =  lLs^(x^y\s). 
This  Az  “diagonalizes”  the  inductance  matrix. 


-4 


-2 


x  0 
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Figure  2.  Three  1000  (am  Wire  Vector  Potential  Az(x,y) 


3.  Variational  Form  of  the  Field  Problem 

Dimensionless  spatial  variables  are  based  on  a  charac¬ 
teristic  time  tQ  of  the  input  pulse  line,  —  ctQ ,  and  a  skin 

dePth  =  \fff- 

With  Sz  —  t0FJv  action  integrals  advance  all  fields, 

if(<D)  =  J dxdA  [(V<£>)2  —  (<9t<J>)2]  ,  (12) 


(VAZ)2-(<9TAZ)2-AZ(-^)VZ  . 

(13) 

Our  finite  element  (FE)  representation  employs  value 
and  time  derivative  degrees  of  freedom,  e.g.  O  = 
[O.  <9tO].  The  variation 

i+iw,-+1  'i+14»+i+1W,.  •,'4>+,+  1  W(._J  ‘“t  =  0  ,  (14) 


yj(  Az  )  =  J  dxdA 


is  equivalent  to  solving  the  PDE  problem.  The  symmet¬ 
ric  positive  definite  weight  matrices  have  the  general 
form,  /  Wm  =  Tlm  ®  (Gx-\-Gy)  —  T(m  0  V  with  integrals  of 
squared  time  values  T ,  time  derivatives  T/,  spatial  gradi¬ 
ents  G,  and  spatial  values  V  taken  over  each  cell. 


C.  Particle  Methods 


The  original  formulation  [3]  used  freely  drifting  parti¬ 
cles  only  and  focused  on  particle  annihilation  and  cre¬ 
ation  rules  to  track  emerging  features. 

The  novel  particle  solutions  used  here  are  drift  kinetic 
fluid  particles  (DKFP)  and  track  precisely  the  distribu¬ 
tion  function  evolving  under  the  action  of  the  drift  ve¬ 
locity  C  and  acceleration  A  which  are  functions  of  local 
position  x,  C  =  V  +  x*5//i,a  =  A  +  x*a//i. 

These  DKFP  contain  three  common  factors  which  rep¬ 
resent  the  dilation  of  the  initial  size  (due  to  the  shear 
in  the  velocity  and  acceleration),  and  the  (asymmetric) 
movement  of  the  initial  domain  boundaries.  For  a  ID  par¬ 
ticle  class  these  are: 


td  1  t2  a 

^)  =  (1  +  T+2— } 

h±  (0  =  3zh  +  (V  +  d=5)  t  +  —  (A  +  d=  cc)  t^ 


(15) 

(16) 


The  expected  profile  for  number  density  (per  unit 
length,  area,  or  volume)  n(X,t)  =<  N  >  ft  is  then: 


V2Ut  J 
(17) 

with  similar  expressions  for  momentum  and  enthalpy. 
The  fluid  properties,  needed  by  the  electrodynamics  to 
support  the  evaluation  of  a  magnetic  Reynolds  number 
(3tm  —  2£Bt0Vr/ 8 '%)  and  the  source  term  above,  are  then 
easily  projected  into  each  node  or  cell  of  the  field  solu¬ 
tion  grid  with  virtually  no  error. 


n(X,t)  = 


N 


2  hD(t) 


erf 


(MOzjO 

V2Ut 


-  erf 


1.  Source  Boundary  Conditions  at  the  Wire  Core 

The  detailed  accounting  of  material  phase  changes  as 
the  wires  vaporize  and  ionize  is  left  to  future  refinements. 
A  reasonably  simple  and  apparently  accurate  approxima¬ 
tion  is  to  balance  the  force  of  the  wire  core  attraction 
with  a  “rocket”  force  due  to  the  inward  flow  of  coro¬ 
nal  plasma [2] .  Here  a  strict  balance  is  too  specialized,  for 
what  is  needed  is  a  path  to  account  for  wire  loads  that  do 
not  hang  at  the  initial  radius  until  they  virtually  disinte¬ 
grate.  From  earlier  work  [4]  we  know  that  A1  foil  switches 
exhibit  an  areal  mass  loss  rate  that  scales  directly  with 
the  surface  magnetic  energy  density,  viz.  dtm  —  —  , 
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and  a  flow  velocity  determined  by  a  sound  speed^l  or,  as 
seen  for  switch  foils,  an  Alfven  speed  characterized  by 
local  field  or  “private  flux”  near  the  wire.  Such  consider¬ 
ations  lead  to  a  mass  loss  estimate,  for  a  wire  of  length  £ 
and  radius  rwire,  that  scales  like 

m  «  10 K{t/rwire)lfMA][iig/ns]  .  (18) 

Typical  values  for  K  are  [5-25]  p  gem /ergs  and  we 
find  that,  in  keeping  with  the  other  models  cited  here, 
only  a  small  fraction  of  the  wire  current  10-2)  can  be 
viewed  as  directly  participating  in  this  erosion  process. 
The  larger  current  fraction  must  be  concentrated  outside 
the  core  in  the  vapor  and  plasma  emitted  by  the  wire. 


2.  Magnetoplasma  Conductivity 

The  Epperlein-Haines[6]  formulation  for  the  required 
skin  depth  80  in  the  wire  corona  can  be  written 


i/tfj _  =  7l1  =  Vo  «l(z> 


(19) 


with,  a ,  (Z,  cor)  a  rational  function,  and 


%  =[^h\  =  1.147519  10 


-14 


L e2neTe 


ZlnA 

J3!2 

W 


H-  (20) 


As  needed  in  the  Az  source  term,  for  a  time  t0. 


(Lf  =  1.218472  - lO-2 


T3/2 


2 10  In  A  Z  a,  <Z,  cot) 


(21) 


is  the  (dimensionless)  equivalent. 


III.  EARLY  INDUCTANCE 
AND  ENERGY  AT  STAGNATION 


chines  like  “Z”.  For  1,2,  and  3  wire  loads  set  at  0.9 
of  the  return  current  radius,  one  sees  quite  favorable 
inductances  for  diameters  near  or  above  100/im.  The 
larger  masses  shown  in  Table  1  above  100/im  clearly 
would  not  be  viable,  but  here  hollow  loads  of  large  radius 
and  lower  mass  should  be  examined. 

With  a  load  mass  mi  in  /igm  per  wire,  inductance  L 
in  nH,  driver  voltage  V  in  MV,  and  pinch  dimensions 
of  length  f|  and  radius  r()  in  cm,  the  available  wire  core 
energy  at  stagnation  scales  as, 

KimPt 2  =  522.4  (r0  (m,i)  l>2  V/L)  [kJ] .  (22) 

One  expects,  Kimp  3  =  \/4j3Kimp  2  for  the  three  wires. 

We  have  (i)  modified  the  wire  dynamic  model  (WDM) 
to  the  calculation  of  kinetic  energy  transfer  to  such  low 
wire  number  loads  on  both  Double  Eagle  and  "Z"  class 
drivers,  and  (ii)  assessed  the  load  energy  and  upstream 
voltages. 

For  DE  any  initial  mounting  of  two  or  three  50  to  200 
gm  Ti  wire  loads  at  radii  in  excess  of  0.6  of  the  return 
current  radius  will  easily  preclude  early  voltages  near  the 
insulator  stack  from  rising  above  one  half  the  open  circuit 
voltage  as  the  load  initiates  and  starts  to  run  in.  The  best 
available  kinetic  energies  are  about  80  kJ  for  two  wire 
loads,  and  about  100  kJ  for  three  wire  loads.  The  (60-80 
pg/ cm)  masses  required  are  typical  for  DE,  the  inductive 
current  “bite”  is  clear. 

For  “Z”  again  initial  mounting  of  similar  wire  loads  at 
radii  in  excess  of  0.6  of  the  return  current  radius  will  also 
keep  early  voltages  in  bounds.  The  best  available  kinetic 
energies  are  about  600  kJ  for  two  wire  loads,  and  about 
900  kJ  for  three  wire  loads.  The  (250-375  /ig/cm)  masses 
used  are  an  excellent  match  in  “Z”. 

In  both  cases  the  result  of  too  much  mass  erosion  is  to 
degrade  the  available  kinetic  energy  at  stagnation. 
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